Time-convolutionless master equation for quantum dots: 
Perturbative expansion to arbitrary order 
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The master equation describing the non-equilibrium dynamics of a quantum dot coupled to metal- 
lic leads is considered. Employing a superoperator approach, we derive an exact time-convolutionless 
master equation for the probabilities of dot states, i.e., a time-convolutionless Pauli master equa- 
tion. The generator of this master equation is derived order by order in the hybridization between 
dot and leads. Although the generator turns out to be closely related to the T-matrix expressions 
for the transition rates, which are plagued by divergences, in the time-convolutionless generator all 
divergences cancel order by order. The time-convolutionless and T-matrix master equations are 
contrasted to the Nakajima-Zwanzig version. The absence of divergences in the Nakajima-Zwanzig 
master equation due to the nonexistence of secular reducible contributions becomes rather trans- 
parent in our approach, which explicitly projects out these contributions. We also show that the 
time-convolutionless generator contains the generator of the Nakajima-Zwanzig master equation in 
the Markov approximation plus corrections, which we make explicit. Furthermore, it is shown that 
the stationary solutions of the time-convolutionless and the Nakajima-Zwanzig master equations are 
identical. However, this identity neither extends to perturbative expansions truncated at finite order 
nor to dynamical solutions. We discuss the conditions under which the Nakajima-Zwanzig-Markov 
master equation nevertheless yields good results. 

PACS numbers: 03.65.Yz, 05.60.Gg, 73.23.Hk, 73.63.-b 
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I. INTRODUCTION 

Electronic transport through small quantum systems, 
such as quantum dots or single molecules, has been inten- 
sively studied in recent years Apart from envisioned 
applications, such devices address fundamental questions 
of non-equilibrium quantum statistics. Quantum dots 
coupled to electronic leads under a bias voltage generi- 
cally relax towards a stationary state. Unless the number 
of relevant degrees of freedom of the quantum dot is very 
small, the relaxational dynamics is complex, including 
broadly distributed time scales and damped oscillatory 
behavior. The stationary state that is eventually ap- 
proached typically depends on the physical parameters 
in a complicated way and can in particular be very dif- 
ferent from the equilibrium state of the isolated dot. 

The descriptions of transport through quantum dots 
or molecules far from equilibrium have so far followed 
three broad approaches. In the first, the focus is on an 
electron tunneling through the device. Its dynamics is 
described by a non-equilibrium Green function (NEGF) . 
The current through the dot can be expressed in terms of 
the local NEGF on the dot, which contains selfenergies 
due to the tunneling or hybridization between dot and 
leads. 4 This hybridization, which is described by a bilin- 
ear component -f/hyb of the Hamiltonian, in principle can 
be incorporated exactly. On the other hand, interactions 
with other electrons, with vibrational modes, or with lo- 
cal spins, which all are particularly important for small 
dots or single molecules, require approximations^— 

The second approach revolves around the non- 
equilibrium Keldysh generating function ! 16 i It is most 
naturally expressed as a functional integral and, with 



suitable source terms, contains the full information on 
the system. This formulation is particularly suitable for 
numerical calculations. When errors due to Trotter dis- 
cretization and a cutoff time for the memory kernel are 
properly controlled, the results are numerically exact. 

The third approach focuses on the dynamics of the 
small system. An equation of motion for the reduced den- 
sity operator in the Fock sub-space of the small system 
is derived by integrating out the lead degrees of freedom. 
The result is a master equation (ME)M^Sr— If the small 
system is sufficiently simple, the interactions within this 
system can be treated exactly. However, integrating out 
the lead states naturally leads to a perturbative series in 
the hybridization -f/hyb- 

Master equations can be either non-local or local in 
time. A non-local ME, for example of Nakajima-Zwanzig 
(NZ) type,2L2£ contains a memory kernel, which relates 
the rate of change of the reduced density operator at 
a time t to the reduced density operator at all previ- 
ous times t' < t. On the other hand, a local ("time- 
convolutionless," TCL) MEi22r— expresses the rate of 
change of the reduced density operator at time t in terms 
of the reduced density operator at time t only. 

If one has a practical method for generating all terms 
in the perturbation series for the transition rates or mem- 
ory kernel in orders of -f/hyb, one can hope to resum the 
series or at least a subseries. This idea has been very 
fruitful for many-particle physics, from the Dyson equa- 
tion to the theory of the Kondo effect. For the non-local 
ME of NZ-type, Schoeller, Schon, and Konig have de- 
veloped a real-time diagrammatic scheme that generates 
all termst ^ 20 ' 25 For a large class of systems including 
a quite general coupling Hamiltonian -ffhyb, Schoeller 25 
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has presented a compact superoperator formulation in 
Laplace space. This formulation is particularly suitable 
for a non-equilibrium renormalization-group approach, 
which in principle includes all orders in iJ hyb i2£ 

Apart from the NZ ME, the T-matrix approach from 
time-dependent perturbation theory has been used to 
calculate the transition rates in the ME. 22 i 32 — It has 
the advantage of being relatively straightforward but 
is known to produce divergences beyond second or- 
der in -ffhyb, the nature of which has recently been 
clarified . 23 ' 26 ' 38 The superoperator derivation of the T- 
matrix ME will make their origin transparent. 

The TCL ME has the obvious advantage of being an 
exact ME describing the full dynamics that is neverthe- 
less local in time. However, so far a method for generat- 
ing all terms in the perturbation series for the TCL ME 
has been lacking, which has limited its usefulness. 

The main purpose of the present paper is to derive an 
iterative scheme for constructing all orders in the per- 
turbative expansion of the generator of the TCL ME. 
The results are valid for the exact ME describing the full 
dynamics. Only at the end we will discuss the implica- 
tions for the stationary state. Furthermore, a surprising 
connection between the TCL generator and the T-matrix 
transition rates is uncovered. This connection introduces 
the divergences of the T-matrix rates into the expansion 
terms of the TCL generator. We will show that these 
divergences cancel order by order. In the present pa- 
per, we concentrate on master equations for the diagonal 
components of the reduced density matrix, i.e., for the 
probabilities. We will call these the Pauli master equa- 
tions or rate equations. 

In the remainder of this paper, the theoretical devel- 
opment is presented in Sec. [TTJ After a brief review of the 
superoperator formalism and the TCL ME, we derive the 
Pauli version thereof, i.e., the TCL rate equations. Then 
we derive the T-matrix formula for the transition rates 
within the same formalism and exhibit the relation be- 
tween the rates derived within the two approaches. After 
that, we organize the perturbative expansion of the rates 
derived within both approaches in a way that makes all 
divergent terms explicit and then show that all the di- 
vergences cancel in the TCL rates. Finally, additional 
insight is gained by a discussion of the Nakajima-Zwanzig 
ME and of the stationary state. The results are summa- 
rized in Sec. IIII1 where we also draw a number of conclu- 
sions. Several proofs are relegated to appendices. 



II. THEORY 

A. Superoperators and the TCL master equation 

Since we will make extensive use of the superoperator 
formalism, we briefly review the superoperator deriva- 
tion of the TCL MEM£2r2±22. To make contact with the 
T-matrix approach and to allow for the analysis of di- 
vergences, we consider a time-dependent hybridization 



between dot and leads. The Hamiltonian reads 

H(t)=H + H hyh e r > t , (1) 

where r\ is small and positive. Hq = i/dot + -ffieads de- 
scribes the decoupled dot and leads. For convenience 
we assume that the eigenstates \m) of i/dot are non- 
degenerated As usual, the leads are represented by non- 
interacting Fermi seas. 

The density operator p of the full system satisfies the 
von Neumann equation 

p=-i[H(t),p]=:-i£(t)p, (2) 

where we have defined the Liouvillian C. The resulting 
unitary time evolution of p can be expressed as 

p(t) = T<_ exp (-i jf dt' C(t')^j p(t ), (3) 

where T<_ is the time-ordering directive. 

Projection superoperators V and Q are defined by 

Vp(t) := [tr lea d S p(*)]<»A eads (4) 

and Q := 1 — V. Here, p° eads describes the leads in gen- 
erally separate equilibrium — each lead is in equilibrium 
at its own chemical potential and possibly temperature. 
We write C(t) = Co + £hybe'' 4 with obvious definitions 
and note the identities 23 



VC = C V, (5) 
VC hyh P = 0. (6) 

We assume that the system was in a product state at 
time to with the leads in equilibrium, i.e., Qp(to) = 0. 



Dropping this assumption would lead to additional terms 
describing the propagation of Qp(to) from time to to t. 
Since we are interested in the case to ~ °o> we do not 
expect the initial state to be relevant. We then obtain 

Vp(t) = PT^exp (-ij dt' [Co + £ h ybe"*']) Vp(t ). 

(7) 

The time-ordered exponential is expanded as 

v P {t) = dt i dt 2 --- dt v 

w= Jto *^*o Jto 

x e- lCo ^-'^hybe^e-^ (tl ~ t2) C hyh e^ t2 ■ ■ ■ 
xC hyh e^e-* Co ^-^Pp(to). (8) 

Equations |(7j) and (|SJ) simply represent the time evolution 
of the density operator, projected with V at time t. 

The TCL ME is derived by splitting the von Neumann 
equation ^ into two parts, 

dVp{t)/dt = -iVC{t)Vp{t) -iVC{t)Qp(t), (9) 
dQp(t)/dt = -iQC(t)V P (t) - iQC{t)Q P {t). (10) 
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The second equation is solved by 

Qp(t) = TVexp (~*Gjf dt'£(t'^j Qp(t Q ) 

-if dt'T^exp(-iQ^ dt" C{t") \ QC(t')V P (t'), 

(11) 

where the first term vanishes under our assumption of 
Qp(t ) — 0. The main idea of the TCL approac h 29 ' 30 is 
to express pit') by propagating the full density operator 
backward in time, 

p(t') = 7% cxp (i J dt" C{t") \ p(t), (12) 

where T_> is the anti-time-ordering directive. Insertion 
into Eq. (flTI) gives 

Qp{t) = -*J t dt'T^exp (-iQ jf dt" C(t")^j Q£(t')V 
x 7% cxp (i J dt" C{t")j [Tp(t) + Qp{t)\. (13) 

Solving for Qp(t) yields 

Qp(t) = [1 - E(f, to)] -1 ^, t ) M*) (14) 

with 

£(t,*o) := dt'T^cxp (-iQ £ dt" C(t")j 

x QC{t')V exp ^ J dt" C(t")^j . (15) 

Inserting Qp(t) into Eq. ([9]) results in an equation of mo- 
tion for V p alone, 

dVp(t)/dt = -iVC{t) [1 - S(i, t )] _1 W)- (16) 

This is the TCL ME for the case Qp(t ) = 0. Using Eqs. 
(HJ) and (JOJ) and noting that £(t, to) contains a projection 
Q to the left, we can rewrite this more specifically as 

dVp(t)/dt = -iVC Q Vp{t) 

- iVC^ 1 [1 - £(*, t )] -1 V P {t) (17) 

with 

S(t,t ) = -iQ f dt' 

x T^cxp (-ij dt" [C + QC hyh e^"Q]^ 

x C^e^VT^exp (i J dt" [C + A iy be" 4 "]^ .(18) 

The first term in Eq. (|17[) describes the unperturbed time 
evolution, while the second stems from the hybridization. 



B. The TCL Pauli master equation 

We here derive an exact TCL ME for the diagonal 
components of the reduced density operator. Since we 
are assuming non-degenerate dot states, this is equiva- 
lent to an equation for the secular part of the reduced 
density operator. 2 - A ME for the diagonal components 
of the density operator, albeit non-local in time, has been 
derived by Zwanzig 2 — and rediscovered, in the context of 
transport, by Leijnse and WegewijSf 2 ^ 

We introduce new projection operators by 



Vp{t) := 



^ |m)(m|tr lcads( o(t)|m)(r 



'PLds, (19) 



where the \m) are the unperturbed dot eigenstates, and 
Q := 1 — P. P evidently projects the density opera- 
tor onto a product form with diagonal reduced density 
operator.— We will call V p the diagonally projected den- 
sity operator. It is easy to show that 



VCo = CoV = 0, 



(20) 



which goes beyond Eq. for V . Since 7/hyb changes the 
electron number in the leads by ±1, we have 



VC hyh V = 0. 



(21) 



We now assume that the system was in a product state at 
time to with the leads in generally separate equilibrium 
and diagonal reduced density operator, Qp(to) = 0. 

Repeating the derivation in Sec. Ill Al with the new pro- 
jections P, Q, we obtain 

dPp(t)/dt = -iVC hyh e vt [1 - t(t,t Q )]- l Vp{t) (22) 

with 



S(t,t ) := -iQ / dt' 

Jto 

x IV cxp dt " [Aj + QC hyh e nt "Q] 

x £ Uyb e"*'p 1% cxp (i £ dt" [C Q + C hyh 



.(23) 



Due to Eq. (|20|). the unperturbed time evolution has 
dropped out of Eq. (f2"2"|) . We can now write 

dPp(t) /dt = S(t,t ) Pp(t) (24) 

with the generator 

S{t, t ) = -iPC hyh e^ [1 - S(t, to)]- 1 ^. (25) 

Equation ([24)) is an exact ME for the diagonally projected 
density operator under the condition Qp(to) = 0. It is 
evidently local in time. Since it only involves the diagonal 
components, we call it the TCL Pauli ME. A Pauli ME 
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in the reduced Fock space of the dot is of course obtained 
by taking the trace over the lead degrees of freedom, 

dpdot/dt = tr leads <S(i, t ) p do t(t) ® Pi eads =: S dot p dot (t). 

(26) 

The reduced generator >Sdot written in the dot eigenbasis 
is the transition-rate matrix. Ensembles of such matrices 
are studied in Ref. |4l| within random matrix theory. 

We have now eliminated the off-diagonal components 
of the reduced density matrix p dot from the equations 
of motion, similar in spirit to Zwanzig's work^S and also 
to Refs. Hj and [26|. We are therefore able to determine 
the dynamics of the probabilities of dot states exclusively 
from the knowledge of these probabilities at a given time. 
This does not mean that we assume the off-diagonal com- 
ponents to be small, which is not generally true. 

The knowledge of the probabilities is sufficient for the 
calculations of dot observables that commute with the 
dot Hamiltonian H do t- To see this, we denote the opera- 
tor for such an observable by A in the Fock space of the 
dot. Then the operator in the Fock space of the whole 
system is A (8 lieads in an obvious notation. The average 
of the observable is 



(A)(t) = T>p(t)A®l 



leads 



trdotPdot(*)^ 



Ti 



(27) 



writing matrix elements of dot states |m), \n) as pf°* = 
(m|pdot|^) etc. If A commutes with the dot Hamiltonian 
we can choose A to be diagonal in the dot eigenbasis 
{|m)} so that 

(m) - £ Ptiit) A mm = TrVp(t) A. (28) 



Thus the knowledge of Tp(t) is sufficient to calculate the 
average. Examples are the charge on the dot, the vibra- 
tional energy of a molecule, or the component of its spin 
parallel to an applied magnetic field, assuming vanishing 
transverse anisotropy. On the other hand, the current 
does not commute with H dot and thus does depend on 
the off-diagonal components^ 2 . However, it is possible to 
reconstruct the full density operator from Vp(t), 

p(t) = Vp(t) + Qp(t) = [I - E(t, h)}- 1 Vp(t), (29) 



compare Eq. ([14 



C. Superoperator derivation of the T-matrix 
formula 

In the following, the T-matrix formula for the transi- 
tion rates is rederived within the superoperator formal- 
ism to allow a direct comparison with the exact TCL ME. 
Moreover, we show that this derivation relies on a single 
straightforward but generally unjustified approximation. 



To start with, note that the exponential time depen- 
dence of the perturbation in Eq. (fTJ) is exactly the case 
considered by Bruus and Flensberg 22 in their derivation 
of the T-matrix formula. The analog of Eqs. (0 and © 
for diagonal projection reads 

Vp(t) =VT^ cxp ^-t jf dt' [C Q + £ h yb e"*'] ) Vp(t ) 



u=0 



dt\ I dt2 ■ ■ ■ dt v C 

to Jto Jto 



hyb 



-'hyb 



Vp(k 



(30) 



where we have used Eq. (l2Tfl) . This is the time evolution of 
the full density operator under the condition Qpito) = 0, 
projected with V at time t. 

Taking the time derivative of Eq. (PU)) , we obtain 



dPp{t)/dt = K(t,t )Vp(to) 



(31) 



with 



1Z(t,t ) :=-tVC hyh e^Y / Hr / dt x / 



ft rtx 

> ' / dt x / dt 2 ■ ■ ■ 

/;,: Jto 



dt a e~ tCo ^-^C^e^e-^ (tl " t2) 



£ hyb e"* 2 ---/: hyb e^P, 



(32) 



where the new summation index is fi = v — 1 and the 
integration variables are now counted by p. The terms 
in this series are of order p + I in £hyb ■ In the limit to — > 
— oo, it is straightforward to evaluate the time integrals 
at non-zero 77, 

K(t, -00) = e unt PChyh [-Co + l^r 1 £hyb 



x [-£ + (v- 2)irj\ 1 £ hy b ■ ■ 
x £ hyb (-C + it])' 1 C hyb P. 



(33) 



The terms are of order v in £hyb- Due to the trace over 
lead states coming from the leftmost V and the equilib- 
rium lead density operator Pi* cads contained in the right- 
most V, this expression contains equilibrium averages of 
products of lead electron creation operators a aka and 
annihilation operators a Q kcr- To obtain a non-zero con- 
tribution, these operators must be paired. Since the hy- 
bridization Hamiltonian iJhyb is linear in these operators, 
only terms of even order v contribute to lZ(t, —00). 

We will now show that the central approximation of 
the T-matrix approach consists of taking 1Z(t,to) to be 
the generator of a Pauli ME that is local in time, 



dPp(t)/dt w n(t,t )Vp(t). 



(34) 



To that end, we show that this ME indeed leads to the 
usual T-matrix formula^ 
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Expressing the Liouvillians in Eq. (|32|) in terms of the 
corresponding Hamiltonians, we obtain 



r 00 ft ft! 

H(t,t )» = -if I ^(-iT / dh dt 2 - 

I p= ■'to 

x H hyh e 7 >\ e~ lHo ( *" tl} [-ffhybe 17 * 1 , e~ tH ° (tl ~ t2) 
x [ff hyb e" t2 , • • • [H hyh e^ , e"** 



x e iH <> (*(•— *o)l . . .l e ii? o t2)j gi-ffo (*— *l) 



dt,, 



(35) 



We now consider unequal initial and final eigenstates, \i) 
and /), respectively, of i?o- Pure initial and final states 
are described by the density operators and |/)(/|, 
respectively. Expanding the nested commutators, except 
for the outermost one, we obtain for the matrix element 
of TZ(t,to) between these pure states 



T fi := (f\{ti(t,t )\i)(i\}\f) 

= -* (/I E H)"*" \ Ml / ^ 



t pt' 



to J to 



dt L 



lit 



dt[ [ 1 di' 2 • • • f " < l^hybe"*, e"^ (t - 4l) 

Jto J t ^ 

x H hyh e^e- iHo ^-^H hyh e^ ■ ■ ■ H hyh e^ 

x H hyh e^e iHo W-^iJhyte^e^ (t "*' l) }|/). (36) 

It is helpful to rewrite this expression as a derivative, 

ft ft\ rtfi—i 



dt 



</iEh)* / dti / * 2 

M=0 "'to Jto 



dt, 



x e 



aHotl H hyh e 71tl e- lHo (tl ~ t2) H h ^ u - ■ ■ ■ 



hybt 



H byh e^e- lHot » \i) 



(37) 



Next, the initial time to is sent to —00 at finite 7/. With 
T\ = t — t\ and = tp-i — t^ for /1 > 1 we obtain 



(/I^THr dn / dr 2 --- dr„ 
1 Jo Jo Jo 



x e 



x e 



- lHoTl H hyh e- lHaT2 H hvh e- lHoT3 ■ ■ ■ e~ iH ^H 



'hybt 



hyb 



-*flb(*— n— TS — — — T>) -T)[ftt— /»Tl— (/i— 1)t 2 — T-J 



dt 



(/IE 



^ Ei — Ef + iptr) 



h 



hyb 



1 



X iJ hyb |i) 



Hi 



1 



hyb 



(38) 



We have used that the \x = term vanishes for |/) ^ 
The fractions are to be understood as inverse ordinary 
operators. The time derivative can now be evaluated, 



r /* = EE 



=1 ^ _ S / _ -E f + ivri) 

1 



X (i|-Hhyb 



x H\ 



H - ir) 
1 



hyb E t -H - iOi - 1)77 

X (/| #hyb 



'Ei-H + i(u-l)ri 



Hhyh\f) 
Hhyb ■ 



1 



H + if] 



H hyh \i). 



(39) 



We notice that the limit 77 — > + can be taken in the 
factors (Ei — Ho ± iftf?) -1 independently from the first 
factor under the sum. In the former, rj > indicates in 
which complex half plane the poles are located. In the 
latter, the limit 77 — > + leads to a 5-function implement- 
ing energy conservation, 



r /i = EE 27r ^- E /) 

fl— 1 z>— 1 

1 



X (t| ^hyb 



hyb 



Ei-Ho- i0+ 
1 

Ei - Ho + i0+ 



■Hhyb I/} 



hyb 



(40) 



Since -Hhyb changes the electron number in the leads by 
±1, Tfi can only be non-zero if \x and v are both even or 
both odd. 

Defining the T-matrix 



T :— 22 -^hyb 



1 



■H 



hyb 



we obtain the well-known result^ 

T fl = 2ir6(E l -E f )\{f\T\i)\ 2 



(41) 



(42) 



Note that we have obtained this result explicitly for the 
exponential time dependence of the hybridization. It was 
not necessary to consider a different time dependence at 
intermediate steps, as in Ref. [H. 

We now use a product basis of unperturbed eigenstates 
m), \n) of the dot and z)), |/)) of the leads. Summing 
over all initial lead states \i)) and final lead states |/)), 
we obtain the T-matrix expression for the transition rate 
from dot state \n) to dot state \m) ^ |n), 

R n ^ m = 2wJ2 W i \((f\(m\T\n)\i))\ 2 S(E n + e l -E m -e f ). 
U 

(43) 

Here, E m (ej) are eigenenergies of dot (lead) states and 
Wi is the equilibrium probability to find the leads in state 
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. The sums over lead states are understood as integrals 
if their spectrum is continuous. 

We have shown that the T-matrix formula (|43|) for 
the transition rates is what one gets if one takes the 
exact time evolution of the density operator, projects 
onto diagonal density operators of product form with the 
leads in equilibrium, and then by hand replaces the pro- 
jected density operator at the initial time, Vp(to), by the 
projected density operator at the present time, Vp(t). 
This confirms the statement made in Ref. [23| that the 
T-matrix approach to transport misinterprets the tran- 
sition rates between dot states \n) at time to — > —oo and 
|m) at time t as transition rates between \n) and |m) 
both at time t. 



D. Relation between TCL Pauli and T-matrix 
generators 

We derive two simple relations between the generators 
S and 1Z. The defining equations (f2"4"|) and (j3"Tj) read 

dfp{t)/dt = S(t,t )Vp(t), 
dfp(t)/dt = K{t,t )Vp{t Q ). 



The first equation is solved by 



Vp(t 1 )=T^expi dt'S(t',t ) )Vp(t 2 ), (44) 



where t\ > ti- Choosing t\ — t and £2 = to arl d taking 
the time derivative we obtain 



^■Pp(t)=S(t,t )T^exp^ dt'S(t',t )) V r ilu}. 



(45) 



Comparison with Eq. (|31j) yields the identity 



1l(t,to)=S(t,t Q )T^exp^J dt'S(t',t )J . (46) 

Conversely, to represent S in terms of 1Z, we integrate 
Eq. (|3"Tj) from time to to t, 

Pp(t) = Pp(t ) + f dt' K(t', t ) Vp{t ). (47) 

Jt 

Comparison with Eq. (|4"4"]l yields 

TV exp dt' S(t', t )j = 1 + dt' TZ(t', to). (48) 



Inserting this equation into Eq. (j46J), we finally obtain 



s(t,t )=n(t,t ) 



1+1 dt'TZ(t',t ) 

to 



(49) 



This remarkable expression allows us to obtain the gener- 
ator of the TCL Pauli ME from the T-matrix generator, 



in principle. This result is potentially useful since we 
have an explicit expression for the transition rates in the 
T-matrix approach in terms of ordinary operators. It 
will also allow us to derive the perturbative expansion of 
S(t,to) in the following. 

The derivation also goes through for the full non- 
diagonal ME. The corresponding expressions can be ob- 
tained by removing the tilde from all symbols. The result 
is equivalent to an identity found by Buzek. 43 



E. Perturbative expansion in the hybridization 

In this subsection we derive expansions of the TCL 
Pauli and T-matrix generators in powers of T/hyb or £hyt>- 
In the following, we send to — > —00 and suppress the ar- 
guments (t, —00). The expansion of the T-matrix gener- 
ator is obtained from Eq. (03]), H = J27=i with 

£(2/») = _ 4 JUVt p Chyh [_£ + ( 2a1 _ l) irl ]-l £ hyb 

x [-£ + (2p - 2)i?7]~ :L £ hyb • • • £hyb (-£0 + iv)' 1 
x C hyh V. (50) 

We have used that all odd orders vanish. 

The TCL generator is obtained from 1Z using Eq. (|49l) . 
The time integral is easily performed, 



h 2m J 



-1 -1 



(51) 



Expanding the inverse and comparing the two sides order 
by order, we obtain 



£(2M) 



Y(-iy y fcpwo^ \ — , 

(52) 

where the second sum is over q + 1 positive integers pi 



adding up to p. We note in passing that Eq. (|52|) can also 
be obtained from the expansion of the TCL generator in 
terms of ordered cumulants, following van Kampen i 39 ' 44 
The first few terms read explicitly 



5(2) = £(2) 



5(4) 
5(6) 



2r] ' 

£(6) _ ^(4) ^_ . 

2r\ 

7^(2) £(2) 



7^(2) 



4?7 



am 



2i] 2r\ 



(53) 
(54) 

(55) 



The first equation shows that in the sequential-tunneling 
approximation the TCL and T-matrix expressions for the 
transition rates agreed 

The problem in exploiting the expansion (|52p is that 
the TZ^ diverge for 77 -> 0+ for all 2p > 4. This is 
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in addition to the explicit divergences due to negative 
powers of 77 in Eq. (|52p . We would much prefer a rep- 
resentation of (S( 2 ^) in terms of expressions that remain 
finite. To obtain one, we first simplify the notation by 
setting t — 0, since in the limit 77 — > + the value of t 
does not matter. We then define 



fc( 2 ^V) -=-if £ 
x [-£ 



hyb + (2/i - 1)177] 1 £ 

(2p-2)ir 1 }- 1 £ hyh ---£ 



hyb 



-hyb 



x (-£ + (2/i' + l)ir 1 y 1 £ hyh V, 



(56) 



where /t > 7/. Note the identity 7v<>> ) = 7v (2 ^. 

Divergences of the type of negative powers of 77 arise 
whenever £q in the inverse superoperators (—£o + inr])~ 1 
can be replaced by zero. These divergences are singled 
out by inserting 1 = V+Q between each pair of £hyb- We 
note that under the assumption of non-degenerate dot 
states, the projection V projects out the secular reducible 
contributions. 2 ® These are thus removed by Q. Since the 
lead-electron creation and annihilation operators must be 
paired between any two V, all expressions with an odd 
number of £hyb superoperators between two V projec- 
tions vanish. Thus at the odd-numbered positions be- 
tween the £h y b, Q = 1 — V does not do anything and Q 
is redundant. This also means that divergences cannot 
arise from the inverse superoperators at odd-numbered 
positions. We therefore only insert 1 = V + Q at the 
even-numbered positions, 

£(2„) = -ifiChyh [-£ + (2/i - l)^]" 1 £hyb (V + Q) 

x [-Co + (2/i - 2)i7?]- 1 £ hyb ■■■(V + Q) £ h yb 
xi-Co+in^Chybfi. (57) 

We denote the regular parts of W^^v) by 

£(2*2,0 . = ^p Chyh [-Co + (2/i - l)^]- 1 £ hyb Q 



x [-Co + (2/i - 2)777] 1 £ hyb • • • Q C hyh 
x [-Co + (2 f i' + l)iT]]- 1 C h y h V, 



(58) 



where a projection Q is inserted at every even-numbered 

position between the £hyb- We also define iticg^ := 

IZieg' ^ ■ The finiteness of Tvlfg' 2 ^ ' for rj — s> + is shown 
in a more general context in appendix [X] Note that 



(2/j,2/x-2) 



= ft.( 2 *2,i-2) and, in particular, izill = it'- 



ll, 

since there is no position to insert Q. This reproduces 
the well-known observation that the second-order rates 
in the T-matrix formalism do not show divergences. 
From Eq. (|57|) we now obtain, using Eq. (|20|) . 



£(3/*) 



-5(2*0) 
/v reg 



£(2,0) 
/Vrcg 

~2v~ 



K 



(4,0) 



477 



-p(4,2) £(2,0) 
/Vreg /vrcg 

4?7 2?7 



-6(2*2,1-2; 

/v reg 



ft 



(2/^-2,2^-4) £(2,(-4,2,i-6) 
' Vrc£ 



ft 



(2,0) 



(2/i - 2)77 (2/i - 4)77 



2n 



(59) 



Since we have inserted V + Q in /i — 1 positions, there 
are 2' 1_1 terms in this sum. In particular, we find 



•75(2,0) 
/v reg ' 



5(4,0) 
/v reg 



5(6,4) 
/v rcg 



£(2,0) 
£(4,2) /x -rcg 

r ° g 277 ' 

-5(2,0) 
' vreg 

277 

r>i 1.2) £(2,0) 



5(4,0) 
£(6.4) /v rcg 

rcg 477 



n 



4t) 2i] 



(60) 
(61) 

(62) 



As an intermediate result, we have thus written the T- 
matrix generator ft order by order in terms of expressions 
that remain finite for rj — ► + and explicit negative pow- 
ers of rj. Since each insertion of V generates a factor of 
I/77, the most strongly diverging term in To 2 ^ scales as 
1/if- 1 . 

Inserting Eq. ([59]) into Eq. ([52]) . we obtain S^ 2 ^ in 



terms of ft 
leading terms read 



(2i/,2iO 
rog 



with < v 1 < v < (i for all p. The 



5(2) 
5(4) 

5(6) 



£(2,0) 

/v Teg 1 



£(4,0) 



£(6,0) 
/v reg 

+ ft (6 

-ft (4 

/v rc( 

-TvL 2 



4) 



5(2,0) 
' v-rcg 

~w 

£(2,0) 
/Vrcg 

277 

5(4,2) £(2,0) 
/vreg /Vrcg 

4?7 2?7 



£(4,2) 
/v rcg 



ft^ 2) 



£(2,0) 
£(2.0) ^ rc S 

r ° S 277 ' 



(63) 
(64) 



£(6,4) 
/v rcg 



ft 



(4,0) 



4?7 



,(2,0) 
"reg 



,(2,0) 



0) 



2?7 277 

£(4,2) £(2,0) 

/Vreg /Vrcg 

4t/ 27/ 



£(2,0) 
_ £(4,0) '^ rc g 

1Cg 277 

£(4,0) 
_ £(2,0) '^rcg 

1Cg 477 

£(2,0) 
, £(2,0) /C icg 
+ 'Vcg ^ 



ft 



(2,0) 
reg 



27? 



(65) 



In this expansion of the exact TCL Pauli generator, all 
singular contributions in the limit 77 — >• + have been 
made explicit. The maximum power is 1/77^ . 

To conclude this section, we illustrate the results by 
considering the terms of fourth order. The corresponding 
term in the T-matrix generator reads 

£(4) = -ifc hyh {-C + 3777)- 1 £ hyb 

X {-C + 2MJ)- 1 £ h yb + 1^ £ hyh f> . (66) 

Let 'RM'' act upon some density operator p. Then 
£hyb (— £a + it])" 1 £hyhVp contains contributions for 
which the second (from the right) superoperator £hyb 
undoes the changes introduced by the first £hyb- Hence, 
£hyb (— £a + ii])^ 1 £hyhPp is an operator with non- 
vanishing diagonal components in the product basis of 
unperturbed eigenstates. But for diagonal components 
we have £ = so that £ in the next super- 
operator to the left, (—-Co + 2irz) ~ 1 , can be replaced by 
zero. We thus obtain a singular contribution proportional 
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to \/2ir}. More formally, we single out the divergent con- 
tributions by introducing 1 = Q + P, 

£W = f Chyh (-Co + Siv)- 1 C hyh (-C Q + 21T])- 1 
x QChy h (-Co + irjy 1 C hyh P 



i PC hyh (-Co + Sir]) 1 £ hy b 



2 i i) 



xW h ,bKo+")) _1 c hyh P 

— P {i) + tz (4) 



(67) 



The divergent part Tt^j*. is identical to Prc'g^ Pilg^/2r], 



according to the definition (|58|) . 

The fourth-order term of the TCL generator con- 
tains a correction term beyond cf. Eq. namely 



/v 



(2) 



277 2 ?7 

= +i VC hyh (-C + ir^y 1 Aiyb 77- 



2,/ 



x 75£ hyb (-/: + ^)- 1 £ hyb 75. 



(68) 



~ (4) 

This looks very similar to the divergent part TZ^'. The 
differences arc the opposite sign and a different prcfactor 
of ir\ in the left-most inverse superoperator. If this factor 
were the same, the correction term would exactly cancel 
the divergent part. As it is, the correction term does 
remove the divergence for 77 — >■ + but leaves a non-zero 
difference behind, 

^dtv - ^ (2) = i^hyb (-A) + ivrH-Co + 3Z77)" 1 
x C hyh Pc hyh (-Co + iri)- l ChybV. (69) 
We will show that this difference indeed remains finite. 



F. Cancelation of divergences 

Our next goal is to show that the divergences described 
by negative powers of 77 all cancel in the limit r\ — > + . 
It is useful to resum the terms in Eq. (f52"j) . 



q=0 (Uo,/il,...,/ig=l 

Inserting Eq. ([5^|. we obtain 



77^(2^1) 77^(2^5) 
2nii] 2\i q r\ 



(70) 



s = y y (-i)"' n^>^o)!i^s — 

(71) 



p— /i ,Mo>A»l!M'nM2,At2>— >Mp 

5(2 M2 ,2^) ^(2^,2^=0) 



2^2?7 



2\i v r) 



where n' + 1 is the number of /4 being zero. The second 
sum is over p+1 pairs (/*,-, fj,^), i = 0,1, ... ,p, with /tj = 



1,2,.. ., /i^ = 0, 1, . . ., and fii > /i-, satisfying either /i- = 
/ii+i or /^ = for any two consecutive pairs. The last 
P'i = ^p mus t equal zero. 

In Eq. (|7ip , 75 represents the explicit order in I/77. How- 
ever, the superoperators Preg' 2fl ^ also depend on 77. To 
find the limit 77 — >• + , we thus have to expand them up 
to the order rf . Their Taylor series in 77 reads 

00 

7^< V) =~i E 

">J|l-l,nl2(l-2,- > m 2(i'+l =0 

x (2/i - l)™ 2 "- 1 (2/i - 2) m2 ^- 2 • • • (2// + 1)™ 2 m'+i 
x [m 2/1 _ 1 ,m 2/1 _ 2 ,...,m 2 ^ +1 ]( 2 ' 1 ' 2 ^), (72) 



where we have defined the notation 
[m 2At _i, m 2At _2, ■ • ■ , ™ 2p ' ■ *l (2 "' 2//) 



lim PAwb 



+ 1 

x [-£ + (2/i - l)^]- 1 -" 12 -- 1 £ hyb Q 

x [-£ + (2/i - 2)i7?]- 1 -"^- 2 £ hyb • • • £ hyb 

x [-Co + (2/i' + lH- 1 -™ 2 ^ £ hyb 75 



(73) 



with Q inserted at all even-numbered positions. In par- 
ticular, Eq. (f72"|) implies that 

lim = - Z [0,0,...,0] (2 ^ 0) . (74) 

77— S-0+ 

It is shown in appendix [5] that the limit 77 — > + in 
Eq. ([73]) converges for all m^-i, m^-i, • • • , TO2p/+i > 0. 
Moreover, we show in appendix [B] that the superoperator 
defined in Eq. ([73"]) does not depend on the values of the 
prefactors of ir\, as long as they are all positive. Thus it 
does not depend on /i and /i' except in so far as 2/i— 2/i'— 1 
is the number of its arguments m v . We therefore drop 
the superscript (2/j, 2/i') from now on. 

Insertion of Eq. ([72]) into Eq. (|71l) leads to an expansion 
of the TCL generator S, 



s = VH) P+1 



E 



(-1)" 



E 



"»0,2>H>-1 TO 0,2,.{, + 1 ">P,2».p-l »»»IM=0 

x 9 ( ~' )Sm 9 (2/io - l)" 1 -™- • • • 
2/ii 2/i 2 ■ • • 2 Hp 

x (2/t + l)™ - 2 ^ 1 • • • (2/ip - l) m p. 2 M P -i . . . 

x [wio^o-i: ■ • ■ : m o,2 M ;,+i] • • • [rrip^p-i, ■ ■ .,m p ,i] 

x ry s ™-P, (75) 

where E m := too, 2mo -i + ■ • • + m Pj i is the sum of all 
m^. The two indices of vrn^ enumerate the factors of 



(2^, 2//) 
reg 



in Eq. (|7ip and the inverse superoperators in 



% 

[m 2ft _i,m 2ft _2, • ■ • , m 2M '+i], respectively. 
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Terms containing positive powers S m — p > of 77 
vanish in the limit 77 — >- + and can thus be disre- 
garded. On the other hand, to obtain a finite limit, 
the prefactors in all terms involving negative powers 
S m — p < must cancel. The cancelations can 
only involve terms with the same superoperator factor 
[mo^o-i, ■ • ■ ,m 0|2(t j + i] • • • [m Pi2flp -i, ■ ■ .,m Pt t]. These 
terms have the same values of p, of the orders 2n^ = 
2j.ii — 2/4, and of all m^. We thus write 



p=0 



00 



no ,ni , . . .,n p — 1 



E H) Em 

m ,2 (10 -i,...,m 02(J j )+1 ,...,mp,2, J p-i,...,mp,i=0 

x /("o, Hi, . . . , n p ; m , 2AI o-i> ■ • ■ 

x [mo,2 W -i, • • - ,"20,2^+1] • • • [mp^np-i, ■ ■ ■ ,m Pt i] 

x rf^ p (76) 

with the prefactors 

/(n ,ni, . . . ,n p ;mo,2 Mo -i, ■ ■ ■ ,m P ,i) 

(-!)«' 



E 



2/i 2 • • • 2/ip 



(2^o - l)™ ' 2 ^-! 



A t o,A'o'A'i>A t i.M2,M2>""A 1 p 

x • ■ • (2fi' + l) m °. 2 ^+! . . . (2fi p - l) m p.^ P -i . . . ; (77) 

where the sum in Eq. (|77j) is now constrained by 2/4, — 
2/4 = 2m being given. With this constraint, the only 
freedom left in the sum is the choice of which (jf. are zero; 
recall that the non-zero /4 equal /Uj+i. The numbers 
/to, Mo, A*i, A*i, A*2, ■ • ■ 1 can be reconstructed from 
the orders no, ni, ■ ■ . , n p and the indices i of the ji! i 
that equal zero. Defining the set 



(78) 



we have p £ Z and n' = |Z| — 1, where \Z\ is the cardi- 
nality of Z . Defining the "non-member function" 



if i S Z 

1 if i & Z 



(79) 



we have tt^ — and n' = p — Yn=a n i ■ Replacing p! { by 
/ii — rii in Eq. (I77|) we obtain 



/(n ,ni, . . . mo^o- 



E 



ns (-1)'' 

2[i\2[i2 ■ ■ ■ 2[i v 



x (2/zo - l) mo ' 2 "0" 1 (2 A i - 2) mo - 2 "o- 2 • ■ • 

x (2u - 2n + l) m o,2™-2„ 0+ i . . . 

x (2^ p -l) m ^ 2 "p- 1 (2^ p -2) m ^ 2 ^- 2 ---l m ^ 1 , (80) 

where fi t — rii + nf n^+i + 7rf 7T l ^_ 1 ni+2 + • • • and fi p = n p , 
and products are understood to equal unity if they do 
not contain any factors. For convenience, we define 



Note that S m = J2i=o M i- 

The evaluation of Eq. (|80|) . which is presented in ap- 
pendix O is rather lengthy but has a remarkably simple 
result: For all S m < p, 

/(n , m,...,n p ; mo^-i, . . . , m Pi i) = 1 (82) 

if and only if 

S ro - p = M + Mi + . . . + M p - p = (83) 

and there does not exist any integer i < p such that 

M + Mi + . . . + Mi - i = 0. (84) 

Otherwise, / = 0. Recall that / is not of interest if 
S m > p. Furthermore, we also show in appendix [C] 
that the condition for non-zero / can only be satisfied 
if m p ,2 Mp -i = • ■ ■ = m p .i = 0. 

Inserting these results into Eq. (|76p. taking the limit 
i] — > + , and renumbering the mi jt/ , we obtain 

00 00 
lim S = -iV(-l) p V 

p—0 no ,ni , . . .,n p — 1 

E 

m , 1 ,...,m 0i 2n -i,---,m p _ ltl ,...,m p ^ li 2n p _ 1 -l 

x 6(M - 1) 6(Mo + Mi - 2) ■ ■ • 
x 9(M + Mi + . . . + M p _i - p) 

X [m ,l, . . . , m ,2n -l] ' ' ' ["V-1,1, ■ ■ ■ , Wp_i,2n p _i-l] 

x[0,...,0], (85) 

where the sum over the m^v is constrained by too,i + 
. . . + irip-i^np-i-i = P an d we have defined 



1 for n < 0, 
W := i 1 forn > 0. 



(86) 



Note that the factor 6(M + Mi + . . . + M p _i - p) is 
redundant. 

The p = contribution in Eq. (1851) does not contain 
any sums over m i v since S m = p = 0. There just re- 
mains a sum over no, the order in -f/hyb, i-e., the p = 
contribution reads — i [0] — i [0, 0, 0] — . . . According to Eq. 

t5(2m) 



Mi := mj )2w -i + mj,2/ij-2 + • • • + m;. 



2^ii-2rii + l- 



(81) 



(|74)l . this equals lim^^ + Z^u^xeg ■ Thus in the expan- 
sion in ffhyb) the expansion term 5( 2 W of the TCL gen- 
erator contains the properly regularized T-matrix term 
7?iog^ plus corrections. Furthermore, all these corrections 
contain [0, . . . , 0], i.e., an expansion term of lZ rcg , as the 
right-most superoperator factor. 

Suppressing the limit directive from now on, we find 
that i S in Eq. ([55)1 is the sum of all terms that can be 
constructed according to the following rules: 

1. Each term is a product of p + 1 = 1,2,... su- 
peroperators of the form [m^i, . . . , rnj,2 nj -i] with 
j = 0, . . . ,p, rij = 1, 2, . . ., and mj^ = 0, 1, . . . 
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2. Defining Mj :— rrij v 



',2n, -i! on ly terms with 



Mq + Mi + . . . + Mj > j for all j < p are allowed. 

3. Only terms with Mq+M\+. . .+M P — p are allowed. 

4. Each term obtains a factor (— l) p . 

We draw a number of conclusions: If an allowed term 
contains a factor [m^i, . . . , rrij^nj-i] then any term with 



this factor replaced by [m' [,.,,,!«', 



i,2n'.-l. 



with 



TOj.i + . . . + rrij^rij-i is also allowed. If 



we denote the sum of all such terms by 



[[ M ]] : =E E [m 1 ,...,m 2 „_ 1 ], (87) 

n— 1 rni + . .. 4-^271-1— M 



we obtain 

oo oo oo 

^ = E(-d p E E 

p=0 A/ = l Mi=max(0,2-Af ) 

oo p-M -Mi-...-Mp_ 2 

E - E 

M 2 =max(0,3-Afo-Mi) Mp_i=max(0,p-A/ -----Afp_ 2 ) 

x [[Mo}}[[M 1 }}---[[M p _ 1 }}[[0]}. (88) 

The last sum is understood to equal zero if the upper 
limit is smaller than the lower one. In order to ob- 
tain an expansion of iS in powers of -f/hyb, we note that 
[mi, . . . , m 2 „-i] is of order H^ h with 2n > 2. Thus [[M]] 
contains contributions of second and higher orders. To 
obtain the expansion term iS^ 2 ^ of order 2p, we thus 
only have to consider terms with p+ 1 < p factors [[Mj]] 
in Eq. (l88l) . The first few terms read 



S( 2 ) = 
= 



-i [0], 

-t [0,0,D]+i[l][0], 
-i [0,0,0,0,0] 
f i [0,0,1] [0] 



(89) 
(90) 



i[0,l,0][0] + i[l,0,0][0] 
i[l][0,0,0]-i[l][l][0]-i[2][0][0]. (91) 



Higher-order terms are easily generated using computer 
algebra. They become increasingly lengthy; contains 
30 terms and S^ 10 ^ already 143. Simplification is possi- 
ble by realizing that some of the terms in S^ 2 ^ contain 
factors of S^ 2 ^ ' with p! < p, as we show now. 

Equation ([55)) is equivalent to the surprising identity 



(92) 



M=0 



The usefulness of this equation rests on the observation 
that [[M]] is of at least second order in i/hyb- Therefore, 
we can express S^ 2 ^ by lower-order terms S^ 2tJl \ p! < p. 
Together with the starting value iS^ 2 ^ = —i [0] = 7?-reg, we 
obtain a recursive scheme for determining S^ 2fl \ 



To prove that Eq. (|92| has the solution given in Eq. 
, we iterate Eq. 



A/ 



iS = [[0]]+ £ [[M ]](-iS) 



M = l 

oo oo 



[[0]]- E E [[M ]][[M 1 ]](- l 5) Mo+Ml - 1 

M = l Mi=0 

oo oo oo 

[[o]] - [[i]] no]] + ]T E E [t M °]] 

M = l Afi=max(0,2-A/ ) M 2 =0 
X [[Mi]] [[M 2 ]] (_ l 5)^+Af 1 +Af 2 -2 



= [[0]] -[[!]] [[0]] + [[M ]][[2-M ]][[0]]T... 

(93) 



M = l 



It is clear how this continues. The terms no longer con- 
taining S are the ones satisfying the conditions in the 
multiple sum in Eq. ([55[) . Thus Eq. (|55|) is a solution 
of Eq. ([§2"|. To show that it is the only solution, i.e., 
that Eq. (|55|) implies Eq. (li?2l , we note that the iteration 
shows that any solution S' of Eq. (|92|) agrees with Eq. 
(|88p order by order in the number p + 1 of superoperator 
factors. Thus we find S' = S to any order p + 1. 



G. Relation to the Nakajima-Zwanzig master 
equation 

The Nakajima-Zwanzig (NZ) ME2L2& and equiva- 
lent formulations are commonly used in the field of 
transport through nanostructures. The real-time dia- 
grammatic technique^— and the suitably generalized 
Wangsness-Bloch-Redfield theory 4 ^— are such equiva- 
lent formulationsi 23 i 26 We again only consider the initial 
condition Qp(to) = 0. The derivation starts in the same 
way as the one of the TCL ME, leading to Eq. ([11]). In- 
serting this equation into Eq. ([9]) and using the identities 
© and dHJ, we obtain the NZ ME^ 

- Vp(t) = -i£ Tp(t) - VC hyh e^ jf dt' 
J dt" (C + QC hyh e 



x T<_ exp 

x QC hyh e^'Vp{t'). 
Expansion in powers of Chyb yields 



r\t" N 



(94) 



u=0 



ft rti 

< I dti / dt 2 ■ ■ ■ 

'to J to 



dt,. 



dt u+ i e 



-iC a (t— ti) 



x QC hyh e^ e~ lCo QC hyh e^ ■ • ■ Q£ hyh e^ 

-iC (t„ 



x e 



' +1, QC hyb e^V P (U +1 



(95) 
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As above, the projections Q at odd-numbered positions 
are redundant, while at even- numbered positions they 
remove divergent reducible contributions . 23 i 26 i 38 

The derivation goes through if we replace V and Q 
by V and Q, respectively. We end up with a Nakajima- 
Zwanzig-Pauli ME for the diagonally projected density 
operator, 



d ~ 



-Vp(t) = -V£ hyh e T > t dt'T^exp 



dt 



i J dt" 



Q£ hyb e"* Vp(t') 



(96) 



The bare time evolution has dropped out because of Eq. 
([2"0)) . It is this ME that is expanded up to fourth order in 
Refs. and[26|. The projections Q now remove only the 
diagonal reducible contributions, not all of them. They 
thus implement the regularization discussed by Roller et 
al£& As in Ref. [26J, the regularization is automatically 
included. Our result shows that it can be formulated 
compactly using suitable projection operators V and Q. 
It has been noted in Ref. |38| and shown explicitly in Ref. 
l26l that the Turek-Matveev schem o 34 i 48 differs from this 
built-in regularization already at fourth order. 

If one is only interested in the stationary solution of 
the ME, Vp{t!) on the right-hand side of Eq. ([96]) can be 
taken to be time- independent. It is then possible to eval- 
uate the time integrals explicitly. The resulting equation 
for the stationary state reads = Q Vp with the genera- 
tor, for to — > —00, 

oc 00 

x [-£ + (2m - I)*??] -1 ChybQ[-C + (2fi - 2)irj\- 1 



£hyb ' ' • Q£hyb (-^0 + "/) 1 Aiyb V. 



(97) 



The redundant projections Q at odd-numbered positions 
have been omitted. Since this is an exact result for the 
stationary state, it should agree with what the exact TCL 
ME predicts. We will return to this point shortly. 

If one is interested in the dynamics, one can still ob- 
tain a local ME from the NZ approach. This requires 
the Markov approximation, which is based on the as- 
sumption that the memory kernel in Eq. (|96[) decays 
rapidly in time. This assumption is often justified since 
relaxation in the leads is rapid but also follows directly 
from the condition of a nearly closed conduction channel, 
I/V < e 2 /h~2- With the Markov approximation, Vpit!) 
is replaced by Vp(t). Taking i — > —00, one obtains the 
approximate "Nakajima-Zwanzig-Markov-Pauli" ME 



dpp(t)/dt = gvp(t) 



(98) 



with the generator Q defined in Eq. (|97l) . 

Comparison of Eq. (|97|) and Eq. (|58|) shows that the 
expansion terms are identical to the properly regularized 



expansion terms of the T-matrix generator (we suppress 
the limit 77 — ^ + ), 



(99) 



Hence, the Nakajima-Zwanzig-Markov-Pauli ME is iden- 
tical, order by order in -f/hyb, to the ME with rates ob- 
tained from the T-matrix approach and regularized by 
dropping secular reducible contributions. Up to fourth 
order, this has been shown by Roller et al£H 

We can now gain additional insight into the failure^ of 
the Turek-Matveev regularization scheme.— The proper 
regularization of the T-matrix expressions can be under- 
stood as omitting all terms in Eq. (|59|) except for the first 
one or, in other words, as omitting all terms contain- 
ing explicit negative powers of r\. The Turek-Matveev 
scheme, applied to the calculation of the fourth-order 
ratesj 3 ^ corresponds to expanding the rates into powers of 
r\ and omitting the diverging part proportional to 1 jr\ and 
then letting 7/ go to zero. The obvious generalization to 
all orders is to omit all negative powers of 77. The two reg- 
ularization procedures thus look quite similar. They are 

not identical, though, since the superoperators Tlrlg' 2 ^ ' 
appearing in the proper expansion (|59l) contain positive 

powers of r\. The positive powers from TZ^g ' 2 ^ ' together 
with the explicit negative powers lead to terms of or- 
der rf , which are retained by the Turek-Matveev scheme 
but are absent in the proper regularization. We reiter- 
ate that both the Nakajima-Zwanzig-Markov-Pauli ME 
and the TCL Pauli ME are automatically regularized — 
for the TCL case this is one of our central results. The 
discussion of the proper regularization scheme is only rel- 
evant if one wants to construct the NZ transition rates 
from the T-matrix expressions. 

The exact TCL ME is not equivalent to the approxi- 
mate Nakajima-Zwanzig-Markov-Pauli ME: As noted in 
the discussion of Eq. (|85p. the p — term in this ex- 
pansion is Y] itrcg^ , which we have now identified as 

the Nakajima-Zwanzig-Markov-Pauli generator Q. Using 
Eqs. ([71)) and ([57| we can also write this generator as 



g = ££g£) = -i[[0]]. 



(100) 



The expansion ([55|) of the TCL generator S contains Q 
as the first term but it is followed by an infinite series of 
additional terms. 



H. The stationary state 

Global conservation of probability implies that a sta- 
tionary solution of any well-formed Pauli ME exists. 
Equation ([M]) then shows that the TCL generator S must 
have a right eigenoperator p s t a t to the eigenvalue zero. 
Due to the V projections in S, this right eigenoperator 
must be of the form 



Pstat 



stat 
Pdot 



1 Pleads 



(101) 
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where pf^t i s diagonal. 

Applying Eq. to p sta t , only the M = term in the 
sum survives and we obtain = [[0]] p s tat, which together 
with Eq. (|100p implies 

Sp S tat=0. (102) 

The reverse is also true: If Q p' stat — [[0]] p stat — then 
Eq. (J88J) shows that 5p^ tat = 0. 

Thus p stat is an exact stationary state if and only if 
/3 s tat is a right eigenoperator of Q to the eigenvalue zero. 
The exact stationary state can thus be obtained from 
the regularized T-matrix or Nakajima-Zwanzig-Markov- 
Pauli generator Q alone, in principle. The formal origin 
of this result is that all corrections to Q in S contain Q 
as the right-most factor, cf. Eq. ([88]). 

There are two caveats, though: (i) The result does 
not apply to approximations obtained by truncating the 
perturbative expansion in i/hyb- It does work trivially 
at second order since = S^ 2 K But already at fourth 
order the TCL Pauli ME for the stationary state reads 

- i([0] + [0,0,0]-[l][0]) pstat = 0, (103) 

whereas the Nakajima-Zwanzig-Markov-Pauli ME is 

- i([0] + [0,0,0]) p^t =0, (104) 

which is not equivalent. 

(ii) The result does not carry over to time-dependent 
solutions. Indeed, if p is any eigenoperator of S to the 
eigenvalue A, Eq. (|92|) gives 

oo 

i\p= ^HA) M [[M]]p. (105) 

M=0 

For A 7^ this does not imply anything for the eigenop- 
erators of Q = — i [[0]]. Conversely, knowing an eigenop- 
erator of Q to a non-zero eigenvalue does not help in find- 
ing an eigenoperator of the TCL generator. For the dy- 
namics, the regularized T-matrix or Nakajima-Zwanzig- 
Markov-Pauli generator Q is not sufficient. 

III. SUMMARY AND CONCLUSIONS 

The dynamics of a quantum dot coupled to electronic 
leads can be described in the master-equation formal- 
ism. To use this formalism beyond the regime of weak 
hybridization between dot and leads, further insight into 
the structure of higher-order terms is required. With this 
motivation, we have derived Pauli master equations (rate 
equations) for the probabilites of dot states, to all orders 
in the hybridization, and both in time-convolutionless 
and time-non-local (Nakajima-Zwanzig) form. Our ap- 
proach uses a projection superoperator V onto product 
states with diagonal reduced density matrix. To fourth 
order, the reduction to the probabilities has been imple- 
mented in Refs. Hi] and [26] by explicitly eliminating the 



off-diagonal components from a Nakajima-Zwanzig-type 
ME. Our approach leads to more compact superoperator 
expressions and is easily generalized to all orders. 

Furtermore, we have presented a superoperator deriva- 
tion of the T-matrix expression for the Pauli ME and 
showed that it fails to take into account the propagation 
of the density operator from the present time t back to an 
initial time to. This answers the question posed in Ref.l23l 
whether it is possible to derive the Pauli master equation 
within the T-matrix formalism instead of using it ad hoc 
to calculate the transition rates. The superoperator for- 
malism has allowed us to establish relationships between 
the TCL Pauli generator S, the NZ generator in the 
Markov approximation (exact for the stationary state), 
Q, and the T-matrix generator 1Z. The off-diagonal com- 
ponents of these generators are the transition rates in 
the respective pictures. Relations between the expansion 
terms of order 2/i, S^ 2 ^ , gK 2 W , and VS 2 ^ , respectively, 
have been given. In particular, the expansion terms 
of the TCL Pauli generator are the sum of the corre- 
sponding terms @( 2 ^ order by order, plus corrections, 
which come from propagating the density operator back- 
ward in time in Eq. (]T2"]) . Only at the second (lowest) 
order the expressions are identical. We have shown that 
both the Nakajima-Zwanzig-Markov-Pauli and the TCL 
Pauli generators converge in the limit r\ — > + , order by 
order. Here, r\ is the rate with which the hybridization is 
switched on. In the NZ case, the absence of divergences 
readily emerges from the superoperator expressions, in 
which the secular reducible terms are explicitly projected 
out, whereas for the TCL Pauli generator it relies on a 
sweeping cancelation of negative powers of r\. 

It is crucial for the derivation that the averages of lead 
operators satisfy Wick's theorem, i.e., that they can be 
decomposed into averages of pairs. Besides reservoirs 
consisting of free fermions as considered here, an anal- 
ogous derivation should be possible for free bosons. 

As is well known, the T-matrix rates diverge for rj — > 
+ . Specifically, the term W 2 ^ diverges as l/?/ 1-1 . The 
divergence noted for the fourth-order term by Averin 49 
thus becomes even stronger at higher orders. We have 
shown that the Nakajima-Zwanzig-Markov-Pauli rates 
g{2n) are idgntica^ order by order, to the T-matrix rates 
with proper regularization. This might lead to an advan- 
tage in practical calculations, as the T-matrix method 
formulated using ordinary operators instead of superop- 
erators is expected to be easier to implement. This reg- 
ularization differs from the one proposed by Turek and 
Matveev^s^ 

As a consistency check, we have shown that the 
stationary state obtained from the Nakajima-Zwanzig- 
Markov-Pauli ME is the exact one, i.e., is identical to 
the stationary solution of the TCL ME. It is quite in- 
teresting how this result comes about: S can be written 
as Q plus corrections that all have S as the right-most 
factor again, cf. Eqs. (]52"]l and (jlOOp . The result does not 
carry over to expansions truncated at finite order or to 
time-dependent solutions, though. 
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Another relation between the expansion terms, Eq. 
|), is also important: all correction terms in S^ 2 ^ 
beyond £?( 2 ^) have factors Q^ 2>1 ) of lower order //' < fj, 
at their right end. Consider the case that all Q^ 2 ^ > for 
p! < fi c are suppressed but (jK 2 ^) is not. Then all correc- 
tions in S^ 2 ^^ beyond the Nakajima-Zwanzig-Markov- 
Pauli term gK 2 ^) are also small. For example, in the 
Coulomb-blockade regime, is suppressed, but is 
not. Then the corrections to the cotunneling rates are 
small in the Coulomb-blockade regime since they involve 
sequential-tunneling rates. On the other hand, deep in 
the sequential-tunneling regime, all fourth-order terms 
are small compared to the sequential-tunneling rates if 
hybridization is weak. However, close to a threshold 
where some crosses over from small to large, for 

example at the Coulomb-blockade threshold, the correc- 
tions can be comparable to the NZ rates. 

Of course, outside of the perturbative regime there is 
no a-priori reason for any term to be small. It is in 
this intermediate coupling regime^ that we expect the 
TCL Pauli ME to all orders to show its power. Since 
the TCL ME describes the dynamics exactly, not just 
the stationary state, and is local in time, it is promis- 
ing for resummation schemes addressing for example the 
dynamical non-equilibrium Kondo effect. 
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Appendix A: Proof of the convergence of certain 
superoperators for 77 — > + 

We want to show that the limit for 77 — » + in Eq. 
(1731) exists and is finite for all non-negative integers 
772>2/i-i) "12^-2) ■ • ■ , W2 M '+i, where jj, > // > 0. It is useful 
to prove a more general statement: For /i = 1,2,... the 
limit for 77 — > + of the superoperator 

VC hyh {-Co + Ca^-lM?)- 1 -™ 2 *- 1 C hyh Q 
X (-£ + C 2M -2«7?)- 1 - m2 "- 2 £ h yb ■ • • 

x £ hyb {-Co + d^)- 1 -" 11 C hyh V, (Al) 

where Q is inserted at all even-numbered posi- 
tions, exists and is finite for all non- negative inte- 
gers 7712^-1, Tniji-%, * ■ * , rni and all positive real numbers 
c 2^-i, C2 M -2, • • • j c%. By a finite limit of a superoperator 
we mean a finite limit of all its matrix elements 



By inserting the completeness relation 

• = Ei^I'I^oi = E<*i-Ij)N)oi (A2) 

ij ij 

in suitable places, the matrix elements of Eq. (|A1[) can be 
expressed in terms of matrix elements of Co, £hyb, and 
Q alone. The matrix elements of Co are 



Tr\i)(j\Co\k)(l\=6 jk 6u (E k - E t ) 



(A3) 



Here, Ek and Ei are eigenenergies of Hq, including dot 
and lead contributions. The former are discrete and, by 
assumption, non-degenerate, whereas the latter have a 
continuous spectrum. The proposition could fail if a zero 
matrix element of Co occured in Eq. (|A1[) . 

At this point it is useful to go over to a single-particle 
description of the leads. As noted in Sec. Ill CI the pro- 
jections V in Eq. (| Al|) introduce equilibrium averages, 
triads ■ • ■ P? C ads> over l ea d-electron creation and annihila- 
tion operators. These averages are non-zero only if all 
lead operators are paired. In the expression (|A1|) . which 
is of order 2/x in £hyb, there are \i such pairs. 

Consider a certain inverse superoperator (—Co + 
c l/ if/ , )~ 1 ~ m " in Eq. (|A1|) . Some of the paired lead opera- 
tors may straddle its position, which is numbered by v. 
For two paired lead operators that are both to the right 
of this position, the superoperator (—£0 + c^ir;) -1- " 1 " 
acts on an operator that is diagonal in the single-electron 
state associated with the paired operators. Its energy 
thus does not appear in the difference E k — Ei in Eq. 
(|A3|) . Consequently, only lead-operator pairs that strad- 
dle the position v contribute to the energy difference. Let 
us denote the number of such pairs by Q u > 0. Then the 
difference E k — Ei has the form E n — E n i + Yli=i ^ e p w >j 
where E n , E n > now denote the energies of the dot many- 
particle eigenstates \n) and |n'), respectively, and the 
Ae p , p — 1,2,..., /j are lead single-electron energies. 
There are /1 such energies, which are independently in- 
tegrated over from —00 to 00. The ordering of the two 
£hyb insertions where the corresponding lead electron is 
created and annihilated determines whether this energy 
enters with a plus or minus sign in the energy differences 
coming from the Co sandwiched between these two £hyb- 
The single-particle energy enters with the same sign in 
all these factors. It is thus possible to absorb all minus 
signs into the definitions of Ae p . 

The integrand in the integrals over Aei, . . . , Ae M as- 
sumes the general form 

F(A £i ,...,A £m ) 

n e ^ - e < + e a 

V=l ^ 1=1 



(A4) 



where the function F contains the remaining dependence 
on the single-electron energies due to Fermi functions and 
possibly energy-dependent densities of states and tunnel- 
ing amplitudes. F is assumed to be a real analytic and 
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bounded function of its arguments. Note that for perfect will remove these, though. 

crystals this does not hold due to the appearance of van We rewrite the expression (|A4|) by introducing two sets 
Hove singularities in the density of states. Any disorder of auxilliary variables x v and A,,, 



2 M -i g ( x » 



F(Ae 1 ,...,Ae M ) J 



dx%--- dx 2fJ ,-i 



2/i 



_ a exp iA^ 



2tt 



2/i- 

"2^" 



2 _ En , +x *+ Culrl y+rn v - ( A5 ) 



The integrand as a function of x„ has a pole of order 1 + m u in the negative half plane. Furthermore, it vanishes 
rapidly for x v — > +ioo [x v — > —zoo) if A„ > (A„ < 0). The only possible exception is the case of m v = and X v = 0, 
which we exclude now and treat separately later. Hence, we can close the integration contour in the upper (lower) 
half plane and obtain 



exp iX v 



dx v 



i=l 



{E n „ - E n , +x u + c v iT 1 ) 1 + m " 





-■^1 {i\ v ) m v exp ( -iX v 



C V VX] 



for X v > 
for A„ < 0. 



(A6) 

Note that for > 1 the case A^ = 0, which is included in both lines, is consistent. On the other hand, for m v = 
this case was excluded. We find that Eq. (|A6I) shows a step discontinuity at A„ — for m u — 0. The result after 
performing the integral over X v does not depend on the value at a single point, though. The expression in Eqs. (|A4[) 
and (IA5I) now becomes 



F(A £l ,...,A6 M ) / dAr 



mi 



m 2 u-i! 



: exp I -i X v E nu - i*^ + ^ Ae Pl/j + c^ir/ 

\ !/=l L i=1 

/•0 _^l+mi AmigCi-jyAi 

F(Aei,...,Ae M ) / dA x ^ exp (-tAi [E ni - E n (\ ) 



,•1+77X2^-1 \" i2 A»-l C2 M -ir;A2 M - 

' ' v 2/i-l ° 



mi! 
exp (-iX 2 f l -i 



dX 



2/i-l 



F„2 



— oo 
At 



exp — j 



e k- 



-+... + A, 



A< 



p i 



(A7) 



r 



where ^~ (z^j~) is the first (last) position for which the 
single-electron energy Ae p appears in the energy denom- 
inators in Eq. (|A4[) . 

Integrating Eq. (|A7|) over all Ae p , we obtain 



_jl+mi^mi e ci7)Ai 

dAi 5- exp (-i'Ai [F ni - E n > ]) 

-oo "ii ! 

_jl+ni2n-l ^™ 2 ' J - 1 gC2,i-l'?A2 tI _i 
^2^-1 



x exp f-iA 2/J _i 



-E 1 ., 



- F„/ 

2jx— 1 



xF(A 1/i -+... + ^,...,A„- + ... + A 4 ) (A8) 
with the Fourier transform 



F( K i,..., Kp ) := / dAe ie - tKlA£l 
x F(Aei,...,Ae M ). 



(A9) 



If there are no lead-operator pairs straddling the posi- 
tion i>, i.e., C, v — 0, the variable X u does not occur in F 
in Eq. (IA8|) . The integral over A„ can then be evaluated 
and is proportional to \c v r\ — i(E, lv — F n / )]~ 1 ~ m ". If in 
addition E nu — E n > vanishes, we obtain a divergence for 
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r] — )• + . But by our assumption of non-degenerate dot 
states, E nv = E n > u implies \n v ) — \n' v ). Thus for this con- 
tribution the dot density matrix is diagonal at position 
v. Because of C, v — we can then insert a projection V 
without changing the result. But there is already a pro- 
jection Q at this (even-numbered) position and we obtain 
VQ = 0. The divergent term is thus removed. On the 
other hand, for E nv ^ E n > there is no divergence. 

It remains to consider the case of at least one lead- 
operator pair straddling position v. Then A„ does occur 
in F. We now consider the properties of the functions 
F and F. The behavior of F at large |Ae p | should not 
affect the transport and we can therefore assume that 
F vanishes sufficiently rapidly and sufficiently smoothly 
for Ae p — > ±oo. We thus assume that all derivatives 
d n F/dAe%, n = 0,1,2,... vanish for Ae p -)• ±oo and 
that all these derivatives are absolutely integrable in Ae p 
over the real axis. These assumptions require the previ- 
ously discussed analyticity property. Under these condi- 
tions we have Hm* -►ioo n p F = for all n = 0,1,2,... 

and all p = 1, . . . , fj,. Thus the Fourier transform F falls 
off faster than any power for k p — ► ±oo for all p. 

It follows that F falls off faster than any power for 
A„ — > — oo. Thus the integral over \ u in Eq. (|A8|) con- 
verges for any m v , E nu — E n > , and rj > 0. It thus con- 
verges pointwise for 77 — ^ + . The convergence is also 
uniform since the integrand in Eq. (|A8I) is bounded in 
absolute value by the integrand in the expression 

/0 %mi f Nm 2(1 -i 

-00 Tf\\ ! J_ 00 m 2M -i! 
x |F(A l , r +... + A I ,+ ,...,A !/ - + ... + A I/ +)| (A10) 

and this integral converges. 

In summary, all terms generated by taking the relevant 
matrix elements of Eq. (jAll) and introducing the com- 
pleteness relation (|A2[) remain finite for r\ — > + . Since 
the number of these terms is finite, the whole quantity 
remains finite. The convergence is uniform. 

Appendix B: Proof of the identity of certain 
superoperators 

To prove the cancelation of divergences in Sec. 
Ill Fl we also need to show that the superoperators 
[m 2 ^_i,m 2jU -2, • • • ,ro 2M '+i] (2 ' 1 ' 2//) defined in Eq. (|73j! do 
not depend on the values of the prefactors of in, as long 
as these are all positive. Therefore, we now prove the 
following statement: In the limit n — > + , the superop- 
erator in Eq. (jAll) is independent of 02^-1, c 2/J _ 2 , • ■ • , c\ 
for all non-negative integers to 2j u-i, 777, 2m - 2 , ■ ■ ■ , mi and 
all positive real numbers c 2jtt _i, c 2/1 _2, • • • , c\. 

As shown in appendix this limit is finite. The 
derivative of (|A1|) with respect to c„ is 

-i(l + m v )n V £ hy b (-Co + Ca^-iir?)" 1 -" 12 *- 1 

x £ h yb Q ■ ■ ■ (-Co + c - vin)- 2 - m » ■ ■ ■ C hyh V. (Bl) 



The derivative and the limit n — > + commute because 
(i) the expression (|A1|) is diffcrentiable with respect to 
c v for all n > 0, (ii) it converges pointwise for n — > + as 
shown in appendix [XJ and (iii) its derivative with respect 
to c v converges uniformly for n + (this is shown by 
a trivial modification of the proof in appendix [X] noting 
that the factor n is bounded by unity for < 77 < 1). 

In Eq. (|Bip . the part V ■ ■ ■ V has a finite limit for rj — > 
+ , as shown in appendixE] Including the extra factor of 
n, the limit vanishes. Consequently, the expression (|Aip 
is a constant function of C2 /i -i, c 2(tt _2, ■ ■ ■ , c\. 

Appendix C: Evaluation of prefactors in the 
superoperator expansion 

In this appendix we evaluate the functions / defined 
in Eq. (|80p . which appear as prefactors in the expansion 
of the TCL generator in powers of 77. We first consider 
the case p = 0, which is more easily done for the original 
expression in Eq. (|77p . This expression does not contain 
any sum since n p = 77,0 is fixed, n! equals zero because of 
fi' p = jj! Q = 0, and we obtain 

/(no;m ,2no-ii ■ ■ • .^0,1) 

= (2/i - l)" l °' 2 "o-i(2n - 2) mo ' 2 "o- 2 ■ • ■ l mo >\ (CI) 

Since we are only interested in the case S m = Mo < p = 
0, the only possibility is Mo = and thus mo, v = for 
all v, giving f(n ; 0, . . . , 0) = 1. 

For p > 1 we evaluate Eq. (|80p by iteration. We first 
perform the sum over ttq ■ The term under the sum is a 
polynomial in It^o^i of order Mo- The zero-order term 
in this polynomial vanishes when the sum is performed 
due to the factor (—1)"° . In particular, for Mq = 
this is the only term and the whole expression vanishes, 
/(n ,ni,...,n p ;0, ...,0,mi )2jUl -i,...,mp,i) = 0. For 
Mq > 1, in all remaining terms of orders 1, . . . , Mo in 
27T(f/Ji, only the ttq = 1 contribution survives. We thus 
obtain a polynomial of order Mo in 2/ii with the zero- 
order term missing. It is thus possible to cancel a factor 
of 2/ji with the same factor in the denominator. What 
remains is a polynomial in 2/ii of order Mo — 1 > 0. 

Now we combine this polynomial with the fac- 
tors (2/ji - l)" ll ' 2 "i- 1 (2/ii - 2) mi ' 2 ^i- 2 •••(2/ji - 2m + 
X)nn,2 M1 -2 rei +i m Eq (|8Qp, These represent a polynomial 
in 2yiti of order Mi > 0. The product is thus a polyno- 
mial of order Mo + Mi — 1 > 0. Using [i\ = n\ + 7rf /i2, 
we obtain polynomials in 27rf 112 of order Mo + M\ — 1. 
If Mo + Mi - 1 = and p = 1, 7rf = is fixed and we 
obtain a non-zero result. If Mq + Mi — 1 = and p > 2, 
we can perform the sum over irf . But only the factor 

(— l)^ 1 depends on 7rf and / vanishes. 

If Mo + Mi — 1 > 1, we necessarily have S m > 2. 
Then we only have to consider p > 2 and there exists 
a sum over 7rf . As before, the zero-order term in the 
polynomial in 27rf/i2 cancels and the other terms only 
survive for 7rf = 1. The result is a polynomial in 2/i2 of 
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order Mq + Mi — 1 > 1 with the zero-order term missing. 
Canceling a factor 2/i 2 with the denominator, we obtain 
a polynomial in 2/i 2 of order Mq + Mi — 2 > 0, which we 
combine with the following term to give a polynomial of 
order M Q + Mi + M 2 - 2 > 0. Analogously to the above, 
if Mq + Mi + M2 — 2 > and p = 2, we obtain a non-zero 
result, whereas for Mq + Mi + M2 — 2 > and p > 3 
we get / = 0. If A/ + Mi + M 2 - 2 > 1, which requires 
S m > 3, we iterate these steps. 

We obtain / = if there exists an integer i < p with 

M + Mi + . . . + Mi - i = 0. (C2) 

We obtain / ^ if this condition is not satisfied and 

M + Mj + . . . + M p -p= E m - p = 0. (C3) 

This implies that M > 1, M + M x - 1 > 1, Mq + Mi - 
1 + M2 — 1 > 1, etc. and thus 

M + Mi + . . . + Mi - i > 1 (C4) 

for all i < p. Finally, if M +Mi+. . . + M p -p = E m -p < 
there must exist an i < p such that condition (|C2[) is 
satisfied and we obtain / = 0. 

We draw some conclusions for the case of non-zero / 
with S m = p. Since M + M x - 1 + . . . + M p _i - 1 > 1 
and M + Mi - 1 + . . . + Mp - 1 = 0, we find M p = 0. 
This implies that M + Ml - 1 + . . .+M p _i -1 = 1. Since 
further M + Mi - 1 + . . . + M p _ 2 - 1 > 1, we conclude 
that M p _i < 1. By iteration we find that Mj < p — i. 

The next goal is to find the non-zero values of / for 
all cases with S m = p. For p = we have found 
/(no; 0, . . . , 0) = 1. For p > 1 we already know that 



"%>,2/i -1 = Tip,2/j -2 = • • • = m-p,i = is required for a 
non-zero result. Equation (|80[) then assumes the form 

/(n , n>i, . . . , n p ; 7710,2^0-1, ■ • ■ > m p-i,2 Mp -i-2ri p _i+i, 

nL'o 1 (-1)^ 



0,...,0) = (-!)» 



E 



2ui2u 2 • • • 2u p 

x (2« - l) mo - 2 ™- 1 ( 2 A*o ~ 2) m °- 2 ™- 2 • • ■ 

x (2« - 2n + l) m »- 2 "o-2»o+i . . . 

x (2u p _! - l) mp - 1 ' 2 ^-i- 1 (2/i p _ 1 - 2) m »- 1 - 2 »p-i- 2 ■ ■ ■ 

x 1 m P -l,2 t , p _ 1 -2„ p _ 1 + l_ (Qg) 

The factors following the fraction contain exactly p fac- 
tors of the form 2(Xi — v = 2(?ij + nfrii+i + irf irf +1 rii + 2 + 
. . .) — 1/ with i £ {0, . . . ,p — 1} and v G {1, . . . , 2n,; — 1}, 
where for > 2 we count 77ij„ factors. We rewrite this 
product as rifc=o (^M-ik ~ v k)i where we assume, without 
loss of generality, < io < ii < . . . < i p -i < p — 1 and 
Vk < if ife = ik' and fc < fc'. Then the condition 
Mq + Mi + . . . + Mi — i > 1 for i < p implies ik < k for 
all k. Thus we have 

/(n , m,...,n p ; m ,2 M -i> • • • > °) 

1 p— 1 ttp -1 ^ ^ 

=(-i) p e n^r- nl r 

^,...,^=0=0 ll*=i^ 
=: / p (n , m, . . . , n p ; io, *x, ■ • ■ , ip-i; ^0, ^1, • ■ • , ^-1), 

(C6) 



where the subscript in f p refers to the number of factors 
2/j. ik — Vk in the numerator. 

By adding and subtracting a constant, we can write 
for any j £ {0, . . . ,p — 1} and any real number c, 

f p (n , ni,..., n p ; i , h, ■ ■ ■ , i p -\] vq, v u . . . , i/ p _i) 

= fp( n o, • • • , ti p ; *o, • • • , *p-i! 

i/ , . . • , fj-l, c, fj+l , • • • , V p -i) 

- i v 3 ~ c ) /p-i("o, • • .,n P ;«o, • ■ -^j-ii^+ii • ■ - 

i/ , ■ • ■ ,Vj-i,v j+ i, . . (C7) 

The second term on the right-hand side contains / p _i, 
which has E TO = p— 1 factors 2fii k — Vk in the numerator. 
But we have shown above that for E TO < p the term / 
vanishes. Thus only the first term remains and we find 
that f = f P does not depend on v.j for any j. Thus we 
can replace Vj by 2(n ij + n ij+ i + . . . + rij) (recall that 
ij < j for all j) without changing the value of /. We 
obtain 

f p (n ,ni, ...,n p ; i , h,. . . ,i p -i;v , Ui, . . . , v p -i) 

= (-i) p E IK- 1 )^ 

^,7rf,...,^_ 1 =0 i=0 



m=i2Mi 



(C8) 



The factor for k — in the numerator contains ik = ia = 
and thus reads 2(no+irQ fii—no) = 2ttq Hi. In the factor 
for k — 1 we have to distinguish the two cases ii = 0,1. 
For ii = 0, the corresponding factor in the numerator 
reads 2(no + 7Tjf /ii — no — n i) — ^{^o n i + 7r 6^ 7r f A*2 — 
This factor is multiplied by ttq from the k — factor. 
Since (71-f ) 2 = 7rf, we can drop the ttq in the k = 1 factor 
and write it as 2(ni -I- 7rf u 2 — «i) = 27rf u 2 . If instead 
«i = 1, the fc = 1 factor reads 2(ni+7rf /i 2 — ni) = 27rf /i 2 . 
We thus obtain the same result in both cases. 
For larger fc, the factor in the numerator reads 

2 (n ik + nf k n ik+ i + ... + Tif^+i ■ ■ ■ ?rf n k+ i 



n ik+ i - ... - n fe ). 



(C9) 

7rf_ 1 from the 

factors for j < fc, we can drop all 7r^ with j < k and 



Since this factor is multiplied by 7^ 

we can drop all 
obtain simply 2ir% jik+i- We finally find 



f p (n , ni, . . . ,n p ; i ,ii, . . . ,i p -V, v , v\, . . . , v p -i) 
1 ttp-1 o z 

=(-d p e n(- i ) 7rf n n° P / fc+1 

= (-l) p (-l) p = 1. (C10) 



We have shown that the coefficients / for E m = p vanish 
if condition (IC2[) is satisfied and equal unity otherwise. 
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